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We present a systematic analysis of all the contributions at the leading log order to the branching 
ratio of the inclusive radiative decay B — » Xs + 7 in the framework of supersymmetry without 
i?-parity. The relevant set of four-quark operators involved in QCD running are extended from 
6 (within SM and MSSM) to 24, with also many new contributions to the Wilson coefficients of 
(chromo)magnetic penguins for either chiral structure. We present complete analytical results here 
without any a priori assumptions on the form of ii-parity violation. Mass eigenstate expressions 
are given, hence the results are free from the commonly adopted mass-insertion approximation. In 
the numerical analysis, we focus here only on the influence of the trilinear X'ijj^ couplings and report 
on the possibility of a few orders of magnitude improvement for the bounds on a few combinations 
of the A' couplings. Our study shows that the Wilson coefficients of the current-current operators 
due to f?-parity violation dominate over the direct contributions to the penguins. However, the 
inter-play of various contributions is complicated due to the QCD corrections which we elaborate 
here. 
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I. INTRODUCTION 

With the successful testing of the gauge sector during the LEP era, the standard model (SM) of particle physics 
has certainly established itself as an essentially correct theory at and below the GeV scale. However, it still leaves 
many basic issues to be investigated. Even if one forgets about the hierarchy problem, and the need for neutrino 
masses, a cursory look at the error bars of about 20 % on the various flavor parameters 0, and their wide range of 
apparently arbitrary numerical values suggests that the flavor sector is still much of a misery. Within flavor physics, 
it is the phenomena of the flavor changing neutral currents (FCNCs) that could be the window of the physics beyond 
SM. Forbidden at the tree level due to the unitarity of V^ckm, FCNCs are loop induced where exotic virtualities (from 
supersymmetry or otherwise) could pop up and hence completely change the predictions for the decay rate. This 
could prove to be complementary to the direct search of the exotic particles by hinting toward a particular class of 
models or by ruling out large regions in the parameter space of beyond SM models, and thus, guiding the accelerator 
searches. 

Among the FCNCs, the process B ^ Xg+j is particularly attractive, both theoretically as well as experimentally. 
From the point of view of theory, though the calculation is very complicated and was a major challenge due to the 
ambiguities following 'scheme dependence' of the decay rate right at the leading log (LL) J^,^], it was worth an effort. 
Unlike many hadronic processes, it can largely be freed from the scale and scheme dependence within the framework 
of renormalization group improved perturbation theory and heavy mass expansion together with the assumption of 
quark- hadron duality. All this is possible because mi, >> Aqcd and hence, the inclusive decay B — > + 7 is very 
well approximated by the corresponding partonic level transition b s + 7 where the non-perturbative pathogens 
play only sub-leading role being suppressed by at least two powers of mt (the corrections are less than 10%). The 
next to leading log (NLL) result for the branching fraction in the case of SM is given as 0, Q : 

Br [B ^Xs+j (E^ > 1.6Gey)]g^,j = (3.57 ± 0.30) x 10"'^ . (1) 

Here, E^ is the energy cut on the photon spectrum to get rid of the background photons. The above result implies 
that this particular 'rare' decay is actually not so rare. This is because, unlike other rare decays, where the rate is 
(X G^aggj^,, here it is oc G^aqEo- Also due to the heavy top in the loop, GIM suppression is not much effective. 
The high rate has already been quite well measured by CLEG 6J, BELLE JJ and ALEPH 0. The results of 
different experiments are consistent with each other. A weighted average of the available experimental measurement 
is problematic, because the model dependence errors (and also the systematic errors) are correlated and differ within 
the various measurements. An analysis taking into account the correlations, leads to the following world average Q 

Br[B X,+-f {E^ > 1.6GeU)]^^p = (3.34 ± 0.38) x 10"^ . (2) 

Within la it matches the SM predictions '3,'^ . So it is clear that there is not much room for new physics contributions^ . 
Either the new physics contributions are negligible, or there are large cancellations among various contributions. Thus, 
a study of new physics contributions can help to reduce the viable regions in parameter space of the relevant model. 
This channel is indeed very well studied in the most popular theory of beyond SM physics, namely the minimal 
supersymmetric standard model (MSSM) 0, HTl . There is already a vast literature on the topic and we are 
not adding to that here. The imposed baryon number (B) and lepton number (L) conservation (in the form of R- 
parity) within MSSM, is totally ad hoc. It is an overkill of a problem of proton decay which can be easily rescued by 
alternatives like baryon parity. It also forbids the neutrino masses that would otherwise arise naturally. In fact, super- 
symmetrizing the SM with its sacred gauge symmetries and renormalizability, one naturally lands into the generic 
supersymmetric standard model or the supersymmetry (SUSY) without i?-parity jl3||. However, due to the presence 
of a large number of free parameters, any phenomenological study of i?-parity violation (RPV) models is a daunting 
task. On the other hand, one may consider the origin and suppression of these parameters as due to a spontaneously 
broken anomalous Abelian family symmetry which is invoked to understand the fermion mass hierarchies [T^ . In 
ref.|l5j such an attempt is made to understand the pattern of ii!-parity violation consistent with phenomenology at 
the electroweak scale. However, in the absence of a compelling model of the kind, it is imperative that all possible 
phenomenological constraints on these, o priori, free parameters be studied. Such complexity necessitates a careful 
examination of the exact definition of the RPV couplings, and the choice of an optimal parametrization of the 
latter, or rather the full model Lagrangian, that can simplify the structure of mass matrices to enable an analytical 
diagonalization in a suitable approximation. The single vev parametrization (SVP), first explicitly advocated in 
is such a parametrization. In this paper, we investigate the RPV contributions from the generic model to the process 



^ See the note added in proof at the end. 
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& ^ s + 7 within the framework of SVP flGj . We give analytical formulae that include all possible contributions 
without a priori assumption on the form of R-parity violation. In the same spirit as an earlier study on /i — )■ e + 7 
|l8j |. we give mass eigenstate expressions, hence free from the commonly adopted mass-insertion approximation. In 
the numerical analysis, we focus in this paper only on the influence of the trilinear A^^j, couplings and report on the 
possibility of a few orders of magnitude improvement for the bounds on a few combination of A^^^. parameters. Our 
formulation also has the advantage that the physical meaning of the A'-couplings would not change at all, whether 
the numerical values of any of the other RPV parameters are vanishingly small or otherwise. 

There have been some studies on the process within the general framework of i?-parity violation. More systematic 
analysis are exemplified by refs. [T^l20l |. Ref. 19], fails to consider the additional 18 four-quark operators which, in fact, 
give the dominant contribution in most of the cases. As we shall demonstrate, the interplay of various contributions 
along with the QCD corrections, is very complicated and hence various contributions could add up or cancel depending 
on a particular case. The more recent work of ref. '20] has considered a complete operator basis. However, we find 
their formula for Wilson coefficient incomplete, and they do not report on the possibility of a few orders of magnitude 
improvement on the bounds for certain combinations of RPV couplings, as we present here. We must mention here 
that around the same time as ref. authors of ref. 21] published an analysis of the CP-asymmetries in the radiative 
B-decays, which involved a similar calculation machinery with the extended operator basis. To the best of our 
knowledge all the available studies on the topic in the literature, in fact are incomplete or incorrect to some extent. 
This is mostly a problem of incomplete consideration for the RPV couplings of the bilinear type and the mass mixings 
the latter produce in the sectors of fermions and scalars. We present here a complete treatment of the full operator 
basis together with the corresponding Wilson coefhcients under the truly generic and consistent formulation of SUSY 
without i?-parity. In fact, the existence of interesting contributions from combinations of bilinear and trilinear RPV 
couplings in the process has been noted in the earlier studies of the model under the SVP formulation [22, HI] . We 
present the full analytical result here. However, due to other complications involved, we will postpone any numerical 
study on the aspect to a later publication. Focusing our attention here on the trilinear RPV parameters, we perform 
our analysis at the leading log order. Such an approximation entails uncertainties of about 25 %. However, owing to 
sensitive dependence of the result on the large number of input parameters, we consider it worth-while to examine 
the relevant parameter space and obtain order of magnitude bounds before performing any precision analysis. 

The paper is organized as follows: In section II we briefly summarize the main features of SVP parametrization and 
also set our notation. Section HI deals with the effective Hamiltonian formulation, the extended operator basis and 
the Wilson coefficients. Section IV discusses the decay rate calculation in relation to the anomalous dimension matrix 
and the renormalization group running of the full set of 28 operators. We discuss our major numerical results in 
section V, and conclude the paper after. Some details involving the full set of 28 operators and the 28 x 28 anomalous 
dimension matrix is left to an appendix. The numerical results discussed here are restricted to those obtained from 
combinations trilinear couplings along. Contributions from combinations of a trilinear and a bilinear parameter, the 
latter going into the loop diagrams through RPV mass mixings, form another class of novel results essentially not 
studied before. We present that latter part, from our parallel study, in an independent publication |24|. 



II. 



THE SINGLE VEV PARAMETRIZATION 



The most general renormalizable superpotential for the generic supersymmetric SM (without i?-parity) can be 
written as 



2 ^apkLa^^pE^ 



\K,kUfDpt 



(3) 



where (a, h) are SU{2) indices, {i,j, k) are the usual family (flavor) indices, and {a, (3) are the extended flavor indices 
going from to 3. In the limit where \jk, ^ijk^ ^ijk ^^'^ vanish, one recovers the expression for the i?-parity 

preserving case, with Lq identifled as Hd- Without i?-parity imposed, the latter is not a priori distinguishable from the 
Li's. Note that A is antisymmetric in the flrst two indices, as required by the SU{2) product rules, as shown explicitly 



here with £12 



1. Similarly, A" is antisymmetric in the last two indices, from SU{3)c- i?-parity is exactly an 



ad hoc symmetry put in to make Lq, stand out from the other L^'s as the candidate for H^- It is deflned in terms of 
baryon number, lepton number, and spin as, explicitly, TZ = (—1)3^+^+25'^ The consequence is that the accidental 
symmetries of baryon number and lepton number in the SM are preserved, at the expense of making particles and 
super-particles having a categorically different quantum number, i?-parity. As mentioned above, _R-parity hence kills 
the dangerous proton decay but also forbids neutrino masses within the model. 

After the supersymmetrization of the SM, however, some of the superflelds lose the exact identities they have 
in relation to the physical particles. The latter have to be mass eigenstates, which have to be worked out from the 
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Lagrangian of the model. Assuming electroweak symmetry breaking, we have now five (color-singlet) charged fermions, 
for example. There are also 1+4 VEV's admitted, together with a SUSY breaking gaugino mass. If one writes down 
naively the (tree-level) mass matrix, the result is extremely complicated with all the /x„ and Xapk couplings involved, 
from which the only definite experimental data are the three physical lepton masses as the light eigenvalues, and 
the overall magnitude of the electroweak symmetry breaking VEV's. The task of analyzing the model seems to be 
formidable. 

Doing phenomenological studies without specifying a choice of flavor bases is, however, ambiguous. It is like doing 
SM quark physics with 18 complex Yukawa couplings, instead of the 10 real physical parameters. As far as the SM 
itself is concerned, the extra 26 real parameters are simply redundant, and attempts to relate the full 36 parameters 
to experimental data will be futile. In the case at hand, the choice of an optimal parametrization mainly concerns the 
4 La flavors. In the single vev parametrizationp^ (SVP), flavor bases are chosen such that : 1) among the Lq's, only 
Lo, bears a VEV, i.e. (L) = 0; 2) y]^:{= Xo^k) = ^ diag{mi, m^, ma}; 3) y^^,{= Xq^^ = ~\jOk) = ^ diag{md, m^, mt}; 



4) /life = -^Vj^KM diag{m„, mc, mt}, where Vg = \/2 (Lq) and u„ = \/2{Hu). Thus, the parametrization singles out the 

Lq superfield as the one containing the Higgs. As a result, it gives the complete RPV effects on the tree-level mass 
matrices of all the states (scalars and fermions) the simplest structure. The latter is a strong technical advantage. 
The soft SUSY breaking part of the Lagrangian can be written as 
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soft 



= Q^ml Q + U'^fn^ U + &ml D + DmlL + &ml E 



ml 



\Hu\ 



-^BB + ^WW 



2 ^ijk^i ^'j^k 



2 ^ijkUi b'^ b'j: 



h.c. 



(4) 



where we have separated the i?-parity conserving A-terms from the RPV ones (recall Hd = ^o)- Note that L^7h?L, 
unlike the other soft mass terms, is given by a 4 x 4 matrix. Explicitly, m^^^ corresponds to m^^ of the MSSM case 
while rhl^Js give RPV mass mixings. Going from here, it is straight forward to obtain the squark and slepton masses. 
No RPV A-term enters the mass matrices for instance. Full matrices have been worked out, together with various 
diagonalizing matrix elements from the perturbative approximations for the matrices within the experimentally viable 
setting of small RPV couplings. Readers are referred to ref. JJJ for further details. 



III. THE EFFECTIVE HAMILTONIAN 



A. The basic Strategy 



The partonic transition b ^ s -f- 7 is described by the magnetic penguin diagram (see FigQ. QCD correc- 
tions are obtained by attaching a gluon line to the quarks lines. Such a correction forms a power series in 
aQCD {Miv)\n{M^/Tnl) ^ 0.7 which is too large to be treated in the perturbation framework. This can be effi- 
ciently handled within the framework of operator product expansion (OPE) and the renormalization group (RG) 
improved perturbation theory p5| . One writes down the effective Hamiltonian (Tics) as a series of local operators 
QiifJ,) multiplied by the corresponding Wilson coefficients Ci{iJL). Both, the operators and the Wilson coefficients are 
function of the scale /i, however, the effective Hamiltonian is expected to be scale independent. Thus, we can write 

Heff = ^ a(M)Q,(M) ■ (5) 

i 

The choice of scale /i is arbitrary, and separates high scale physics (> /i) and the low scale (< fx). It is typically taken 
to be the mass scale of the decaying particle (in this case the 5-quark) . There are three major steps in the evaluation 
of the decay rate. They are: 

• Matching at Ci are first evaluated perturbatively at the scale at which one decouples the heavy modes 
(say at /i = M„,). This is done by imposing the equality of the effective, and the underlying theory Green 
functions, at external momenta that are much smaller than the masses of the decoupled particles. 

• Resummation: Since the theory contains two disparate scales mb and My, the large logarithms need to be 
resumed using RG improved perturbation theory. The RG equation for the Ci are given as: 



(6) 
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FIG. 1: Magnetic penguin diagram with Higgs in the loop. Here the chirality fiip is on the external 6-quark. Note that the 
photon can be attached to any electrically charged particle inside the loop but not outside the loop. Diagram with a photon 
outside the loop does not contribute to the magnetic transition. There are similar diagrams with other SUSY particles in the 
loop. For the chromo-magnetic penguin one attaches a gluon, instead of a photon to the quarks in the loop. 




'^R 

FIG. 2: A magnetic penguin digram but with the chirality flip inside the loop. 



Here ^ji is the anomalous dimension matrix (ADM). It is obtained from the ultra-violet (UV) divergences in 
the theory. The derivation of ADM is the most difhcult part of the calculation. Solving the RG equation one 
gets the Ci at the scale mfc. Note that in our case the resummation is performed at the LL order. 

• Matrix elements: Having obtained the Ci{mb) one finds the hadronic matrix elements of the relevant operators 
(Qi)- Typically this involves some non-perturbative method; but in case of inclusive B decays, non-perturbative 
corrections are sub-leading and one can approximate the decay by a corresponding partonic level transition. 

B. The Operators 

Having outlined the strategy, let us list the army of operators beginning with the SM 0. So far as one neglects 
the QCD corrections, the relevant operator is the magnetic penguin. However, once the penguin is dressed with the 
gluonic lines, it mixes with other current-current and QCD-penguin operator. They are given as; 

Ql.2 = {sLaYbLa,!}) (cLff Jt_L CL(3.a) , (7) 

Q3,4 = {sLaYbLa,!}) ^ 7/^ 9L»/3a ) , (8) 

i—u,c,d,s,h 

Qbfi = {sLaYbLocp) ^ {qRifilt,qRil3,a) , (9) 
i—u,c,d,s,b 

Qr = TT^ SLa O-piy &flQ F^'" , (10) 

28 = -^mtSLo.<J^..bn0t-"^'G^^''' . (11) 

Here, we have used Fierz re-ordering to write the first two operators, a and /? are the colour indices and are the 
generators of SU{3). Other notations are self explanatory. Operators Q3 to Qe are the QCD penguins, and Qi,s 
are the magnetic and the chromomagnetic penguins respectively involving up-type quark and W-boson. We would 
like to remark here that for a novice there is some chance of confusion regarding the choice of operator basis. The 
above mentioned operator basis is arrived at by considering the external states on-shell and hence these cannot be 
transformed into each other using equations of motion |26l |. In an off-shell calculation one must consider additional 
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operators ^. Even for the on-shell basis there exist one more basis choice in the hterature introduced by the authors 

of m 3. 

Since the magnetic penguins are chirahty flipping operators, they are necessarily proportional to the fermionic mass 
term. In principle one can also write down the following two operators obtained by flipping the chirality of Qy.g- 

27 = -^m,SRo,(J^,,hLc.F^"' , (12) 
Qs - Y|^m,Sfl„a^,6L^i'^"'^G'^'^''. (13) 

However, these are suppressed by a factor of mslvrn, and hence generally neglected in the SM analysis. They do 
give significant contribution in case of SUSY without i?-parity under the discussion in this paper though. MSSM 
contributions are ver y w ell studied in the literature. We shall make only brief remarks over here and refer the reader to 
the literature [TolflTllT^i for a detailed discussion. In MSSM, over and above the usual, W-t loop contributions, one can 
have charged Higgs-up-squarks, chargino-up-squarks, gluino-down squarks and neutralino down-squarks contributions. 
Since there is left-right scalar mixing in the squark sector, there is a possibility of a mass-insertion inside the loop 
(see for instance Fig(21 as example with a mass- insertion inside the loop) . Since these are mass-terms of the already 
decoupled heavy modes, they are included in the corresponding Wilson coefficients that encode the dynamics of heavy 
modes. So in the MSSM, there are new contributions to the Wilson coefficient of Qy, the operator basis is usually 
assumed to be the same as in SM. However, this is not quite true. It has been demonstrated in that there could 
be 100 extra scalar, vector and tensor-type four quark operators, formally higher order in strong coupling (compared 
to the SM operators). The 40 vector type operators are generated by two gluino box diagrams and gluino-gluon 
penguin diagrams. However, the vector four-quark operators do not mix with the SM operators at one-loop and hence 
contribute only at the next to leading log. The 60 scalar and tensor type operators are generated by the two gluino 
box diagram. These could in-principle mix with the Qifi- However, this mixing turns out to be sub-dominant. The 
scalar and tensor type four-quark operators mix only with the dimension-6 penguins (corresponding to the chirality 
flip on the external quark) and hence their contribution to the amplitude is suppressed by mi,/mg relative to that of 
dimension-5 ones ^SOj. So we shall neglect all this extra 100 operators. 

Giving up the restriction of i?-parity adds to the list of four quark operators. From an inspection of the super- 
potential it is straightforward to write down the relevant operators. It is clear that there cannot be any contribution 
from the A couplings. One gets 12 new four quark operators from the combinations of two A'-couplings and 6 new 
four-quark operators from the combinations of two A"-couplings. With so many extra operators of different origin, 
a notation that is intuitively suggestive, compact, unambiguous and at the same time offers convenient expression 
of summation in formulas is virtually impossible. We make a simple choice for the notation as follows: Keeping the 
notation for SM operators as a yard-stick, we divide the operators according to their chiralities (the advantage of such 
a classification would soon become clear). Those operators that lead to 6^ transition are appended to the list of 
SM operators whereas the operators leading to £t — s- transition are denoted with a tilde (i.e., in the form Qi) and 
form a new list. Obviously the above two classes of operators are related hy L R flip. This choice is not without 
potential ambiguity. For instance, the notation Qg, Qig typically denotes semi-leptonic operators, but since we are 
not dealing with semi-leptonic process in this paper, there should not be any ambiguity. With these comments we 
write down the additional operators obtained from RPV. 



A'A' case: 



Qg-ii = {sLaj'^bL/i) {qRff^^tqRa) ; q = d,s,b, (14) 

Q9-13 = {sRal^bRf}) {qLf3lt,qLa) ; q^d,s,b,u,c, (15) 

QSA = {sRal^bRa,!}) X! {qRipY' qRiP,a) , (16) 
i—u.c.d,s,b 

Qbfi = {sRal^bRcn) ^ {qLipl^ qLil3,a) ■ (17) 



^ For instance see the ref.Q for an off-shell basis. Ref. l27l gives a careful treatment about the issues regarding equations of motion and 

effective Lagrangian, often neglected in other works. 
^ The operators are written in such a manner that one does not encounter ill-defined trace of the product of odd number of gamma 

matrices including a 75. This allows them to use the fully anti-commuting 75 in dimensional regularization. 



7 



A" A" case: 



Sl4,15 = (sRal^ bRa,p) (u Rf3 ^ p. U R/3 ^a) , 
Ql6,17 = {sRal'^bRa.p) {dR(3 dRff.a) ■ 

Thus the total number of operators becomes 28 which mix upon QCD renormahzation. 



(18) 
(19) 
(20) 



C. Wilson Coefficients 

The most important operator for the decay 5 — > s + 7 is of course the magnetic penguin Q7. In order to obtain the 
corresponding C7 we first need to define the relevant interaction Lagrangian. At this point we would like to emphasize 
that since we are working in the most general setting of RPV, there is an important distinction in the way we write 
down the interaction Lagrangian and the Wilson coeSicients as compared to ref. 20|. Our analytical formulas contain 
all possible contributions from RPV including the bilinear RPV couplings. In ref. j2^] the treatment of bilinear (RPV) 
parameters is slightly erroneous. They work in a basis where bilinears (/x^'s) have been rotated away. This does not 
get rid of all the RPV effects beyond those described by the trilinear RPV couplings. RPV mass mixings are still to 
be found among the fermions and the scalars f^-ct such mass mixing effects through the "sneutrino" VEVs, 

render the definition of the flavor basis of the charged leptons and the down-sector quarks in the formulation of the 
Lagrangian ambiguous |l3j . The problem is totally avoided in our formulation. Our formulae are given in terms of 
exact mass eigenstates, instead of the often used mass insertion approximation. Because of the RPV mass mixings 
our formulae for the Wilson coefficients do not factor into the MSSM contributions and RPV contributions as treated 
in ref. . With these remarks we list here the interaction Lagrangian. 



Gluino-quark-squark vertex 



= .g,*(.g)$t(J„ 



h.c. 



where, 



GZl 



2 ■^(i+3)m 



(21) 



(22) 



with V Ad^V = diagjAl"} while is 6 x 6 being the down-squark mass-matrix. Here i,j,k = 1 to 3 (same 
below), m = 1 to 6 for the scalar (squark) mass eigenstates. 

Chargino(charged-lepton)-quark-squark vertices 



Ct 



1-75 , 1+75 



*(d,) +h.c. , 



(23) 



where. 



IK 
92 



(p+3)m 



92 



A' 



— ~~ ^CKM '-^2ri J^hm ^ ''cKM '^(q+2)n J^hm 



92 



(24) 



Here, and throughout the paper, the m index counts the scalar mass eigenstates and the n index that of the fermions. 
Matrix diagonalizes the 6x6 up-squark mass matrix i.e. I?"^A^^I?" = diag{A^u}, ^^'-^ m = 1 to 6 for the 
up-squarks. Similarly V^McU ~ diag{A^c} for the n = 1 to 5 charged fermions (charginos and charged leptons). 

Neutralino(neutrino)-quark-squark vertices 



1-75 



-AC!. 



1+75 



(25) 
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where, 

K'm. = -\/2 {tan^(ai - T3/)Xi; + Ta/J*^;} Vt - 2?(1+3)™ 
A'* 

AC:. = \/2 tan^ Oi Xi„ V^.^^^^ - jQ^vi, ^ Jl21 p^c*^ , (27) 

9 2 .92 

where T3/ = —5. The neutral fermions indexed by n have the 7x7 (neutrahno-neutrino) mass matrix to be 
diagonalized by X. 

Charged Higgs(slepton)-quark-quark vertices 



1-75 _^ 1+75 



*(ci,) +h.c. 



(28) 



where, 



9i 



X. 



/-•R' W<*i T/m* -T-)( I T/Tlfe* -T-)i 



52 



52 



CKM ■'^(j+2)m • 



(29) 



I?' being the diagonalizing matrix for the 8x8 charged-slepton charged-Higgs mass matrix. 



Neutral scalar(sneutrino)-quark-quark vertices 



/:^"=52M/(d„)$t(o 



*(di) +h.c. 



(30) 



where 



^1 



^/2 



2to 



7mJ 



pin 



V2 



p+2)m 



in 



(p+7)mJ 



92 vz 



A', 



92 V2 



q+2)r 



(g+7)mJ 



(31) 



being the diagonalizing matrix for the 10 x 10 neutral scalar mass matrix. 
Next we present all the Wilson coefficients at the scale The most important is of course C7. It can be 

decomposed into various contributions as follows. 



Ct = C^ + + + Cf +Cf 



(32) 



The terms from left to right are the SM, gluino, chargino(charged-lepton), neutralino (neutrino), charged 
scalar (sleptons and Higgs) and neutral scalar(sneutrino) contributions. All except the last one are common to MSSM 
too. Note, however, that in accordance with the above formulation, we are putting together in C7 and C7 contribu- 
tions from the light fermion mass eigenstates, namely the charged leptons and neutrinos. Below we list the individual 
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contributions for a decay dj —^di + J. {j = 3 and « = 2 for 6^5 + 7) . 

,2 \ /^2 



= 

C| = 
Cf = 

cf = 

Cf = 

cf = 



3,g| ™1 




(33) 
(34) 

(35) 
(36) 



(37) 
(38) 



Each term includes, implicitly, summation over the n fermion and m scalar mass eigenstates, except that the unphysical 
Goldstone modes are not to be included. The latter contributions are incorporated into the gauge invariant C^ result. 
Contributions to C7 except for the case of C^ are obtained by replacement L ^ R. The C^ term above is given in 
terms of explicit couplings, rather than the effective coupling vertices. In accordance with the notation of other terms, 
the vertices involved in the C^ term above should be two 'L '-couplings the 'i?' counterparts of which vanish. Any 
term above of the LL type, i.e. with two 'L'-couplings has, like C^ term, a chirality flip from the 6-quark external 
line. The rest of the terms, of the RL type, have chirality flip inside the loop as illustrated also by the explicit fermion 
mass ratio factor. Contributions from terms of RR type to C7 coefficient needs chirality ffip from the s-quark external 
line and arc neglected here, as also in most of the literature. One can similarly map out the details for the chirality 
structure for the C7 coefficient. We note, however, that the latter receives extra contributions from yet another class 
of diagrams with A-^j.-couplings and the chirality flip on the 6-quark. Such contributions are missing in the coefficient 
C7, except for the ones suppressed with the chirality flip from the s-quark. 



Cf = 



4M^ ^n/ji^nfei^ft,+3,m-^fc+3,m 



-QdFi 



QuF2 



M? 



-QuFi 



Mi 



M? 



(39) 



Similarly one can obtain the expressions for the Cg and Cg by introducing the colour factors at the relevant places. 
The loop-functions F1-4 are the Inami-Lim functions, given as: 



Fi{x) 
F2{x) 
Fsix) 
Fi{x) 



I2{l-xf 
1 

12 (1 - xY 
1 



{2 + ^x-&x'^ + + &x\nx) , 
{l-Qx + dix"^ + 2x^ -Qx^ \nx) 



2(1 -x)3 
1 

2(l-a;)3 



(-3 -I- 4 a; -a; -2 Ina;) , 
(l-x^ + 2x \nx) . 



(40) 
(41) 
(42) 

(43) 
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CKM *CKM ' 



Below we write down the non-zero effective Wilson coefficients for the current-current operators. 

C2 

1 



C9 
Cio 
Cii 

C9 



V2 
1 

1 

1 



1 



1 



8M 



-A 



2~^V23 



8M 



2 ^i32 



(44) 
(45) 
(46) 
(47) 
(48) 
(49) 
(50) 
(51) 



where, 



~ 8Af^ •^il2^j"13-'^2+i,m-'^2+j",m ) 



Ci4 

C'le 



1 



\'* T/"S 

2 ^il3 '^CKMi/s-'^2+i,m-'^: 



2,m •> 



^y.2/fc2?2,™^2:™^CKM^C 



1 



8Af, 
1 



\'* T/cs „, 7-)i 7-)/* 

,f2 ^j23 ^CKMi/s-'^2+i,m-'^2,m : 



8rn 



2 ^li2^lj?,'-^i+j,'m'-^3+i,m 



-C15 , 



8mi 



"'^ii2^7"i3-^; 



'jl3 3+j,m 3+i,m 



"t 



— —^17 , 



(52) 



(53) 
(54) 
(55) 



(2?2+.,™ + *^?+vm) (252,^+^25^™)* ■ 



(56) 
(57) 



Here, ys and yt denote the strange and the bottom Yukawas. Again the sum over repeated indices is assumed and 
the unphysical Goldstone modes are to be dropped from the sum. Note that the terms proportional to Yukawas are 
missing in the ref ^ ,20.] . This could lead to interesting contributions with the A' on one vertex and a SM Yukawa on 
the other vertex with a RPV mass mixing effect hidden within the V expression. In the language of mass-insertion 
approximation, these kind of contributions involve a RPV mass- insertion {e.g. from a /ii or a Bi coupling) along the 
scalar propagator. 
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IV. THE ANOMALOUS DIMENSION MATRIX AND RG RUNNING 

The RG evolution from the scale M,.,- to the scale appropriate for the decay dynamics (rrib) is described by eq.lO- 
This requires the knowledge of the anomalous dimension matrix 7^ (ADM). Derivation of 7^ presents many subtleties 
and was a main obstacle in the consistent calculation of the decay rate about a decade ago |3, 1^ ■ Below we mention 
the salient features of the calculation. 

• Since QCD does not know the sign of 75 there is no mixing among the sets of Qi and Qi operators. Thus the full 
basis of 28 operators splits into two invariant sub-spaces of 11 and 17 operators related hy L ^ R inter-change. 

• An operator of a dimensionality n can mix into the operators of dimensionality < n. Thus the current-current 
four quark operators influence the RG evolution of (chromo)magnetic Wilson coefficients, but not the other way 
around. 

• At the one-loop level current-current operators do not mix with the (chromo)magnetic penguins. Thus, one has 
to evaluate this mixing at the two loop-level while still working at the leading log approximation. Because of this, 
the ADM and hence the Wilson coefficients at the LL are found to be regularization scheme-dependent. Such 
a scheme-dependence cancels with the corresponding scheme-dependence of the possible finite one-loop (but 
0{a^)) contributions from certain four-quark operators to the matrix elements of 6 — s- s -f 7 . Such contributions 
vanish in any four dimensional as well as t-Hooft-Veltman regularization scheme but not in naive dimensional 
regularization. Such a situation is taken care of by expressing the ADM and the Wilson coefficients in a scheme 
independent manner 31]. 

• Since one needs to perform a two-loop calculation, one must properly take into account the contributions from 
the so called evanescent operators (that vanish in the limit D = A dimension) while working in the D A 
dimension. 

With these remarks in mind one can formally write down the decay amplitude as: 

11 17 

^ = C7 ( S7 I Q7 1 tree + (S7 I Q7 1 &) tree + XI ( '5'^ I 2» I ^) 1 -ioop + XI ( S7 I S» I ^) 1 -ioop ■ (58) 

i=l 1=1 

In order to obtain a decay rate that is regularization scheme-independent, we must express the decay amplitude A in 
terms of scheme-independent effective coefficients. Decay amplitude A is written as: 

A^Cf {s-f\Q7\b)tree + Cf (s7 | Q7 1 fo)tree ■ (59) 

To derive the effective coefficients one makes the observation that in whatever scheme one chooses to work, the one- 
loop matrix elements of the current-current operators for the decay 6 — > s -|- 7 is always proportional to the tree-level 
matrix elements of the (chromo)magnetic operators. Thus, one can write 

< S7|Qj|6 >i_loop = 2;i < S7IQ7I5 >tree, 

< SgluOn|Qi|6 >i_loop = < SgluOnlQsl^ >tree , (60) 

< S7|2»l& >l-loop = Vi < S'jlQrlb >tree , (61) 

< sgluon|Qi|6 >i_ioop = < sgluonlQsl^ >trco • (62) 

In order to find out the proportionality factors yi,Zi etc., one has to evaluate the finite contribution due to the 
insertion of certain four-quark operators into the diagrams with a closed b quark loop (in the NDR scheme) . The only 
operators that could give finite contribution are the ones with the chirality structure (LL)(RR) or (RR)(LL). Within 
SM basis, these are Q5, Qq. With RPV you could also have Q5, Qe and Qn, Qn. The results for {yi}, {zi}, {?/;} and 
{zi} are given as: 



-1 i = Qe, Qii 
otherwise 



1 « - Qs 
otherwise 



(63) 
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With this one can write down the scheme-independent Wilson coefficients as 

11 

1=1 
11 

j=l 

_ _ 17 ^ 



i=l 
17 



i=l 



(64) 



The index i above runs over only the current-current operators. 

The RG evolution of the effective coefficients is determined by the scheme- independent effective ADM. The RG 
equations are given as: 



In fi 



(65) 



Here C^^ are the components of a column vector C'^^ which contains coefficients Cjf as well as C^^ . Since QCD does 
not know the sign of 75, these two sets of coefficients do not mix. The effective ADM at leading log is obtained as: 



-(O)cff 



if 







off 



(66) 



7°*^ represents QCD mixing of 11 operators whose chirality structure is similar to those of SM operators. 7°*^ represents 
the mixing of 17 operators whose chirality structure is obtained by L ^ i? replacement with SM like operators. We 
present the explicit matrices in the appendix 1X1 

In general, the solution of the RGE for the Wilson-coefficients is given by 



C^'^i^i) = V 



c<=^(M^) 



(67) 



J D 



where V diagonalizes (7°*'')^ 

if = V-\-i^^ fV . (68) 
Pq = 23/3 is the one-loop beta-function and "jf is the vector containing the eigenvalues of 7°*^. In our case 



46 



-16 


-16 


-16 


-8 

28 


-8 

28 


-8 

32 


3 


3 


3 



-16 
4 

32 
3 



-16 

4 

2.233 



(69) 



2.233 6.266 6.266 -13.791 -13.791 -6.486 -6.486) . 
With as{Mz) = 0.121 and /i = mf, = 4.2 GeV, the coefficients Cf{mb) and Cf{'m}j) are given as: 

Cf{mb) = -0.351 C^^(Mv^)-h 0.665 ^^^(Mw-)-^ 0.093 ^^^^(Mw)- 0.198 Q^(Mw 
-0.198 Cl^{M„) - 0.178 Cl^iM^) , 



Cf{mb) = QMlCf{M^) + QmbCf{M^) + QmiCf{My,)~Q.l%SCf{My,) 



-0.198 C^o (Mv) 
+0.381 Cll{M^) 



0.178 Clf{M„) + 0.510 Cl^{M„) + 0.510 Clf{M„) 
0.213 Cl^{M„) . 



(70) 
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The numbers multiplying the different Wilson coefhcients are all of the same size, hence no term can be neglected a 
priori. 

Finally, the branching fraction for Br{b — s- s + 7) is expressed through the semi-leptonic decay b — > u\cev so that 
the large bottom mass dependence m^) and uncertainties in CKM elements cancel out. 



Br{h ^ s + 7) = 



r(6^g + 7) ^ 

1 (0 ^ u\ceh'e) 



(71) 



where Brcxpib u\ceve) = 10.5%. Here the presence of RPV couplings offers new decay channels and new con- 
tributions to the semi-leptonic decay rate. In principle, lepton flavor violating couplings demand a summation over 
the three (anti)neutrino mass eigenstates (assuming the ncutralinos are heavy) before squaring the amplitude. In 
considering the A'-coupling contributions, the dominating neutrino mass terms would be at one-loop level. As all the 
neutrino states have close to zero masses, and very small admissible mixings with the gauginos and the Higgsinos, we 
neglect these effects here and stick to effective massless electroweak neutrinos. The summation is then taken simply 
over the three families. Under the same spirit, we neglect here the effects from the lepton flavor violating couplings 
that could change the electroweak character of physical electron. The partial widths in the expression above are given, 
at the LL order with m^/m^ set to zero, as 



r(6^ 57) = 

T{b — > u\c e D) = 

fie) = 
A 

Br 



647r4 



mX 1 



{/(e) + + + iSip + ICip} , 



1927r3 32 
1 - Se^ + 8e^ - - 24e4loge , 

2V2ai^CKM , 

3 



y— 



yrn 



E 

i=l 



van r^d^ T^d -trrm 



e^mn3^2+m,l^2,l 



V™ r - 1 2 



Here e = rric/mi,. The above expression for the semi-leptonic decay rate is the most general one that also includes 
the contribution from bilinear-trilinear combination of _R-parity violating parameters (the term oc he) which has been 
generally missed in the earlier studies. Such contributions have been shown to play a very important role as elaborated 
in detail in our parallel report on 6 — > s -I- 7 and for the case of quark dipole moments in |23| . 



V. NUMERICAL RESULTS ON THE TRILINEAR COUPLINGS 



Our presentation of the numerical results in this paper restricts to the effects of the combinations of trilinear 
couplings only. We present some details of the numerical results on the phenomcnologically interesting A'-couplings. 
The A"-couplings are found not to give interesting results — a feature we will comment on at the end of the section. 
To set the stage for the more detailed discussion on our numerical results, we first recall the salient features of this 
particular decay within the MSSM. Within MSSM, the story is already quite complicated. In the limit of exact SUSY 
there is no contribution to the decay rate due to the cancellations among particles and sparticles [s^. With SUSY 
broken, there could be gluino, chargino, charged Higgs and neutralino contributions depending on what particle is 
exchanged in the loop. Chirality flip can be induced by mass-insertion on the internal or the external fermion line. 
The squark mass matrices are not necessarily flavor diagonal and the flavor violation is introduced as family off- 
diagonal mass insertion in the loop. The couplings of quarks and squarks to the neutral gauge bosons and gauginos 
are flavor diagonal. The more exact treatment, as adopted here, is to use the mass eigenstates for the particles 
running in the loop and hence have the flavor violation effects absorbed into the effective couplings which contain 
the corresponding elements of the squark diagonalizing matrices. With the inclusion of RPV couplings, there are a 
lot more flavor violating couplings, and more new types of admissible diagrams. These are described in our mass 
eigenstate expressions most naturally by simply enriching the effective couplings. 

The significance of the MSSM contribution to the decay rate is highly dependent on the choice of values for the 
background parameters of the model. Over and above the delicate issue of cancellations among various contributions, 
it is also a function of 'how much' and 'where' does one introduces the flavor violation. The approaches in the 
literature can be broadly classified into two categories, (a) The so called constrained MSSM (cMSSM) where one 
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starts with the universal boundary conditions at the high scale and derives the weak-scale spectrum by RG evolution. 
CKM angles are then the sole source of flavor violation 33J. In such a scenario, the dominant contributions are due to 
the chargino and charged Higgs loops. Gluino and neutralino couplings being flavor diagonal, rely on flavor violating 
mass mixing terms in the down squark mass matrix generated from RG evolution. Such mixings within cMSSM 
are typically very small to contribute to the decay rate. Gluino contributions suffer from hard GIM cancellations. 
While neutralino contributions can escape GIM cancellations through Yukawa couplings but cannot compete with 
chargino contribution that are enhanced by heavy top in the loop 11]. (b) Another end of the spectrum is simply 
the unconstrained scenario. Here, the idea is to assume arbitrary structure for the squark mass matrices and derive 
the weak-scale bounds on the off-diagonal flavor violating entries by requiring to reproduce the decay rate within the 
experimental uncertainties 34]. We adopt the philosophy of unconstrained scenario studies. However, we will switch- 
off the i?-parity conserving flavor violations, except that from the CKM matrix to focus on the effect of the RPV parts. 
In fact, we will focus on a pair of A'-couplings (or A"-couplings) at a time while the other flavor violating parameters 
that are otherwise not theoretically constrained will be switched off. Note that the flavor violations originating from 
A'-couplings are essentially of supersymmetric origin though SUSY breaking, and phenomenologically viable, scalar 
and gaugino masses are needed to avoid the intrinsic SUSY cancellations. 

Before our analysis of RPV contributions, we have first checked our code for the case of MSSM in order to reproduce 
the trends and features obtained in the previous works. Here, the matching can be done at a qualitative level only. Most 
of the detailed numerical results in the literature for the case of MSSM are worked out within the cMSSM scenario. 
The low energy spectrum obtained by RG evolution and a particular SUSY breaking boundary condition contains 
correlations among various parameters and there are RG induced flavor violations mixed with that from the CKM 
matrix. With this limitations, our code does reproduce qualitatively, trends and features of various contributions 
obtained by others |0, 0, hj. We give an illustration of our MSSM results in Figs l3l4l and |31 For a scenario 
characterized by squarks around 300 GeV, sleptons around 150 GeV, down- type Higgs and A parameter 300 GeV, 
(up-type Higgs mass and Bg parameter determined from potential minimization condition), gaugino M2 — 200 GeV 
(assuming the relationship Mi = 0.5M2,M3 = 3.5M2) and = -300 GeV for Fig H and 300 GeV for FigEl (note 
that the MSSM parameter ^ is denoted by /Zo in our model notation), we plot branching fraction versus tan/3. As is 
well known, the charged Higgs contribution is of the same sign as the SM whereas the sign of chargino contribution 
depends on the sign of the product Atfi (see Fig. 0] where we have plotted various contributions to the Wilson 
coefficients). In the large tan/3 region, the decay rate is completely dictated by the chargino contributions and hence 
the enhancement or the suppression of the rate actually depends on whether Atii is positive or negative, respectively. 
Within the framework of mSUGRA (see ref.T^), choosing the gaugino masses to be positive, forces At to be negative 
at the weak scale, and hence one requires a positive /i for the rate at large tan/3 to fall within the experimental limits. 
In Fig l3l5l we plot the decay rate for negative and positive /j, respectively with At fixed to be positive. One can see 
in the figures that the decay rate does fall initially with tan/3 and then rises in case of negative fi but for positive 
/i, due to constructive interference between chargino and charged-Higgs the rate rises with tan/3. This can be easily 
compared with the Fig. 2 of ref.j^^j for example. 

While investigating the influence of RPV we kept the spectrum similar to the one selected for the case of MSSM 
with sign of /i negative and tan/3 — 37. Note that in the large tan/3 region QCD corrections drastically suppress the 
rate at the scale mt, (see for instance FiglSJ. This is expected because, in the large tan/3 region, the Wilson coefficient 
C7 is negative due to large negative chargino contribution and at the scale rrib after the QCD correction there is a 
destructive interference between C2 and C7. The choice of spectrum is guided by requiring that in the MSSM limit 
we do obtain branching fraction well within experimental limits. Now the purpose is to investigate the behaviour 
of branching fraction once we turn on the knob of i?-parity violation. As we will see, both, enhancement as well as 
suppression is possible depending on which operator(s) are playing the major role. Both, the direct contribution, 
through the magnetic penguin, and the effect induced by QCD renormalization, through the various four-quark 
operators are significant. 

Due to the large number of RPV parameters (48 just from the superpotential, with more from the soft SUSY 
breaking sector), without further knowledge of their plausible range of values it is impossible to extract any useful 
information on a phenomenological study, with all of them playing a role at the same time. Hence, the sensible thing 
to do at this stage is to focus on a minimal number at a time and study their possible impact. For the case at hand, 
a pair of parameters is taken at a time, while all the others are switched off. We discuss first the A'-couplings. There 
are two possible combinations in which two non-zero A' can give a finite contribution. These are (A) Xi3jXi^2k 

(B) \'j3^tik2- illustrate their effects in more detail below. 

We will compare below the probable bounds we obtained with the similar bounds on the relevant combination of 
A'-couplings in the literature. We have, first, to bring to the readers attention that the existing bounds in the literature 
are typically obtained assuming a sparticle spectrum of around 100 GeV. However, such a low lying spectrum is very 
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FIG. 3: Branching ratio as a function of tan/3 in MSSM with At positive and n negative. The '+' sign is for the QCD corrected 
rate and the 'x' sign for the rate without QCD corrections. Two horizontal lines in this Fig. and all the following figures for 
branching fraction is the experimental uncertainty at la. See text for the values of various parameters. 



dangerous for 6 ^ s + 7 ^ and quite unrealistic at least for the squarks. In consideration of that, and for a better 
comparison with the available 6^5 + 7 result for the i?-conserving contributions (of the MSSM), we stick to a 
slightly heavier spectrum of 300 GeV squarks. Hence, the existing bounds must be rescaled by a factor of three for a 
better comparison with our bounds. Moreover, for some the combinations that contribute to 5 ^ s + 7 there do not 
exist any direct bound. We have then made use of the best bounds available on the individual parameters to obtain 
a bound on such a combination. 



1. Case A: \'i-ij\'h2k 

If the indices i,j,h,k are un-constrained, one can form 81 combinations of two A'-couplings that can in principle 
contribute. However, in the limit there is no flavor off-diagonal mixing in the (s)neutrino mass matrix, the combination 
with i ^ h will not contribute. The case j ^ k requires extra source of flavor off-diagonal squark mixing to contribute; 
hence is also not of interest here. So we shall confine ourselves to the case i — h and j — k and hence study the 
impact of only nine (9) combinations, picking one of them to be non-zero at a given time. Note that so long as 
one sticks to the phenomenologically required suppression in all the flavor off-diagonal mass mixings from SUSY, the 
RPV contributions which require such mixings would be doubly suppressed and hence highly unlikely to be of any 
numerical significance. For the magnetic penguin operators, the contribution is of LL type. 

The above combination will lead to a 6fj — *■ transition with nib mass-insertion on external line. The contribution 



* To be exact, the top squarks contribute directly through the chargino diagram without the CKM suppression and have to be heavier. 
The bottom and strange squarks contribute, if their flavor violating mixings are nonzero. They are hence likely to be required to be 
heavier too. Otherwise, we can actually live with only heavier stops but a light sparticle here. Such a split in sparticle spectrum is 
considered unlikely as it goes in the opposite direction as predictions from most available SUSY breaking models. 
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FIG. 4: Various contributions to Wilson coefficients versus tan/3. SM (line with symbol '+') and charged Higgs (denoted as 
CS-L, line with sign) are of the same sign (positive) whereas the chargino contribution (denoted as CH-L, line with 'x' sign) 
has a strong tan/3 dependence and is of opposite sign (negative) because At^J, < 0. 



can come from chromo(magnetic) penguins (only sneutrino and neutrino loops are possible) or from the current- 
current operators Qg.io.ii- Sneutrino loops being proportional to the inverse of light slepton mass, are dominant 
compared to the neutrino loops which are suppressed by heavier squark mass. The important feature here is that the 
RPV contributions interfere with the SM and the MSSM contributions. Of the possible nine combinations, we will 
discuss three representative combinations with best bounds. The other six combinations follow similar pattern, as 
their contributions have a similar structure. 

In our numerical analysis, we take real and equal values for two A'-couplings at a time. As discussed above and 
illustrated by the analytical formulae, the contributions depend on the complex product. The strategy here then 
simplifies the results to one-parameter scenario for easy presentation, without compromising the physics. Consider 
the infiuence of the combinations |A']^3j.A']^2fcl- F'or fc = 1,2 the physics is identical because the relevant Wilson 
coefficients (Cg^io) have identical structure and so we shall discuss the combination |A']^32Ai22l- FigEl we have 
plotted the branching ratio against (— -^'122 )• Interestingly, the rate falls below the experimental limit (at 

^^'132 = '^'122 — 0.07) with increasing A' value, reaches a minimum and then rises again (the rise is not shown in the 
fig.). To understand this behavior we have plotted all the relevant Wilson coefficients against A' values in Fig[71 It 
can be clearly seen that Cr{Miv) is falling with increasing A' value. This is because of destructive interference between 
the negative i?-parity conserving contribution (sum of W-boson, charged Higgs and Chargino loop) and the positive 
sneutrino loop contribution, as seen in the figure. At the scale rub after QCD running there is again destructive 
interference between the dominant tree level contribution from Cio and C7 [see eq. H70|l ]. because of which Cj^mb) 
further falls with increasing A' value. So because of these destructive interference, the branching ratio is falling below 
the experimental limit and we obtain a bound of |Ai32Ai22l < 4.9 x lO^'^, to be compared with a rescaled existing 
bound of 3 X 10~^. In FiglSJ we have plotted the branching ratio against A333 {— A323). Unlike the above case, 
here the relevant current-current Wilson coefficient Cn has a more important role to play. Not only does the scheme- 
independent effective coefficient influence CT{mi,) through RG running, it also contributes to Cy^ directly. The 
effect cancels with the i?-parity conserving contributions and hence C^^{Mf^) is seen to reduce with increasing A' 
value. Even at the scale mf,, after QCD running there is a destructive interference between the i?-parity conserving 
and i?-parity violating parts and hence Cy^ at mt, further reduces resulting in the branching fraction that falls with the 
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contributions, respectively. The letter 'L' indicates that these contributions go into C7. Similar notation would apply to the 
following figures. The letter 'R'(L) would indicate that the contribution goes into coefficients with(without) tilde. 

RPV coupling. The fall is much stronger compared to the previous case. We obtain a bound of IA333A323I < 1.6 x 10"^ 
to be compared with the rescaled existing bound of (1 — 3) x 10~^ jS^. 



2. Case B: A,%Al 



Such a combination leads to a 5^ — > transition for 771^ mass insertion on the external line. This would mean that 
the h and s quark at the effective vertices would be right handed. It contributes to the Wilson coefhcients of several 
operators. These include the coefficients of the (chromo) magnetic penguins, via charged lepton, charged slepton, 
neutrino and sneutrino loop. Interestingly, in the SM and MSSM, this coefficient does not receive much contribution. 
The above combination can also contribute to Wilson coefficients of the current-current operators, Qg to Q13. These 
current-current contributions are mostly the dominant ones but not always as we will soon see. 

As discussed at length in case (A), we must require i = h and j — k and hence case (B) then effectively consists 
of the combinations ^'*j3^ij2- Consider Fig 1 101 where we have plotted the branching fraction against non- vanishing 
^233 (= ^232). We get a bound IA233A232I < 1.9 x 10^^, to be compared with the rescaled existing bound of 4.4 x 10^^ 
[33. The mild behaviour of branching fraction suggests the existence of cancellations. In Fig^]we have plotted 
branching fraction against A'j^^a (= ^112) ^^'^ obtain a bound of |A']^]^3A2]^2l < 1-4 x 10~^. Although the branching 
fraction behaves more or less similarly in both the above cases, the underlying dynamics is quite different. To better 
appreciate the differences, we have plotted the relevant Wilson coefficients in Fig^]for the combination IA233A232I 
and in Fig^jfor the combination |Ai]^3A']^]^2l- Fig ll2l the horizontal line (circles) is the large i?-parity conserving 
contribution to C7(mh) which does not interfere with the RPV contributions going into C7. The surprising feature 
is the large (though less than i?-parity conserving contribution) positive contribution of C7(M„,). This is surprising 
because it receives dominant negative sign contribution from the charged-slepton (crosses in the figure) whereas the 
positive sign sneutrino contribution (star marks in figure) is sub-dominant. Here again, the requirement of scheme- 
independence plays an important role because Cj^{ii) — Ct{ii) — Cii(/^). It is the large negative sign contribution 
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FIG. 12: Relevant Wilson coefficients as a function of A233 {~ A232). 

from Cii{Mnr) which adds up with the positive sneutrino contribution and dictates the behaviour of Cj{Miv). At 
the scale uib after QCD running, there is some suppression in the Wilson coefficient. Consider the case of the 
combination |A'ii3A'ii2l (see FigEl- It has a pecuUarity of the charged-scalar loop contribution shghtly dominating 
the contribution from current-current Wilson coefficient Cg. This suggests that it is not always true that tree- level 
contribution dominates. Such a dominance can be traced to two effects : (a) In contrast to the tree level contribution, 
a charged slepton loop contribution has contributions from two diagrams corresponding to photon being emitted from 
a charged slepton (cx to loop function F2) and the photon being emitted from an up- type quark (cx -Fi). (b) The values 
of loop functions depend on the mass of the quark in the loop - the lighter the quark less the suppression. Thus, 
provided that there is no CKM suppression, up-quark-charged-slepton loop dominates over the top-quark-charged- 
slepton loop and in fact can also dominate over the tree level contribution. That there is no CKM suppression in both 
the above cases can be easily verified from the formulae of the Wilson coefficients. Also the uniform dominance of 
charged slepton loop over sneutrino loop can be accounted for by the argument (a) above. Since there are two possible 
diagrams for charged slepton loop it brings in a charge factor of (2/3 -I- 1 = 5/3) which is five times bigger than the 
charge factor coming from the down-squark in the sneutrino loop. However, the charged-slepton and sneutrino loop 
contributions always come with opposite sign and hence there are cancellations. 

In Tabled we have listed bounds that we have obtained for the various A'-coupling combinations among the cases 

(A) and (B). As far as one restricts oneself to the situation of two non-zero trilinear A'-couplings, case (A) and case 

(B) are the only combinations that can contribute. With non-zero bilinears as well, or additional non-zero trilinears, 
there could be more combinations that can contribute. For instance one can have a combination of A^g^ A^j^jji gi'^ii^S 
rise to a neutrino or a neutral scalar loop for the penguins. Since the overall process conserves lepton number, the 
above combination obviously needs to be supplemented with two more L-violating couplings. This could be from 
a Majorana mass-insertions for neutrino or sneutrino in the loop. However, such a contribution would be further 
suppressed and hence we skip presenting explicit results of the kind. 

In all the above discussion, we have assumed that all the couplings are real and positive. It would be interesting 
to consider the possibility when there is a relative phase resulting in a relative negative sign among the chosen non- 
vanishing pair of A' couplings. It turns out that such a case does not affect the bounds much. For the case (B) 
with dominant contributions from RPV couplings, it obviously cannot affect the decay rate which is modulus squared 
of Wilson coefficient. However for case (A), a negative or a positive sign can lead to constructive or destructive 
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TABLE I: Bounds on the various products of A' obtained compared with the appropriately rescaled existing bounds |3^ 



Product |A'A' 


Our bound 


Existing bound 


Wilson coeff. 


Case A 








131.121 


4.9 X 10"^ 


2 X 10"^ 


C7, Cg 


132.122 


4.9 X 10"^ 


3 X 10"^ 


C7, Cio 


133.123 


1.6 X 10"=* 


7.2 X lO"-'^ 


C7, Cii 


231.221 


4.9 X 10"^ 


2 X 10"^ 


C7, Cg 


232.222 


4.9 X 10"=^ 


2 X 10"^ 


C7, Cio 


233.223 


1.6 X 10"^ 


2 X 10"^ 


C7, Cii 


331.321 


4.9 X 10"^ 


(1 -3) X 10"^ 


C7, Cg 


332.322 


4.9 X 10"=* 


(1 - 3) X 10"^ 


C7, Cio 


333.323 


1.6 X 10"^ 


(1 - 3) X 10"^ 


C7, Cii 


Case B 








113.112 


1.4 X 10"^ 


1.3 X 10"^ 


C7, Cg 
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1.4 X 10"^ 


1.3 X 10"^ 


C7, Cio 


133.132 


1.9 X 10~^ 


7.2 X lO"'^ 


C7, Cii 


213.212 


1.4 X 10"^ 


1.3 X 10"^ 


C7, Cg 


223.222 


1.4 X 10"^ 


1.3 X 10"^ 


C7, Cio 


233.232 


1.9 X 10~^ 


4.4 X 10"^ 


C7, Cii 


313.312 


1.4 X 10^^ 


1.3 X 10"^ 


C7, Cg 


323.322 


1.4 X 10"^ 


1.3 X 10"^ 


C7, Cio 


333.332 


1.9 X 10"^ 


(1.3-9.2) X 10^^ 


C7, Cii 
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interference between i?-parity conserving and violating contributions. A relative negative sign though changes the 
sign of the slope for the branching fraction it does not alter its magnitude. Thus, the only implication as compared to 
the earlier bound which (for instance) was due to a branching fraction rising above the experimental limit, would now 
change to a bound due to the branching fraction falling below the experimental limit. The change in the magnitude 
of the bound is however nominal. 

Finally, we comment briefly on contributions from A"-couplings. We obtain no useful bounds on the A"-couplings. 
The main reason behind this is that squark spectrum of 300 GeV is heavy enough to suppress the A" contributions. 
However, a lighter spectrum (although not viable due to large MSSM contributions) cannot lead to enhancement 
because of certain cancellations taking place. This provides a nice example of the interesting role played by the QCD 
running and hence worth a brief discussion. The combinations responsible for contributing to the decay rate can be 
written as: ^h*nz^kn2 from the up-squark loop and ^''*h'^^nk2 from the down squark loop, thus contributing to C-j^. 
Depending on the particular value for the indices, there could also be contributions from the current-current operators 
with Wilson coefficients, Ci(= — C2),Ci4(= — Ci5),Ci6(= — C17). It might appear that the relative negative signs 
for the above coefficients leads to cancellations and hence there is no impact on the decay rate. This is not what is 
happening. There are cancellations, but these result rather from the form of C7 at the scale (m^) (see ea. (|70f) . Taking 
into account, the antisymmetry in the last two indices, one can form 18 combinations from the above mentioned general 
form of the A" contributions (12 combinations contributing to the down-squark loop and 9 to the up-squark loop but 3 
three combinations common to both). There are two reasons that kill the contributions. In a few cases, the loop level 
contribution is accompanied by current-current Wilson coefficients pairs C14 and Cig , or Ci and Cie . Although these 
pairs contribute about equal magnitude with identical sign, at the scale mf, after QCD running these are multiplied 
by about equal coefficients but with opposite sign and hence the cancellations. The rest of the combinations get 
contributions from only one among the above pairs and hence cannot suffer from above cancellations. In this case, 
they all require flavor violating mass-insertions for (s)particlcs in the loop. The only exception to above two cases is 
the combination X'l*^X'l^2 fo'^ which we get a bound of 0.5. This bound is very weak when compared to existing bound 
from other sources ( 6.2 X 10-3[33). 

VI. CONCLUSION 

In summary, we present a complete analysis of the decay rate 6 s + 7 at the leading log order for the generic 
supersymmetric SM, or SUSY without i?-parity. Unlike previous studies on the topic, our results are fully generic, 
admitting all possible forms of i?-parity violation without a priori assumption. We use exact mass eigenstates in 
our formulae, hence free from the otherwise commonly used mass-insertion approximation. In case one prefers, our 
formulation does provide perturbative approximations through which the explicit dependence on all RPV parameters 
can be extracted [13 ■ We consider the analytical results useful for any detailed study on the model in relation to the 
radiative B decay. 

In the numerical implementation of the analytical formulas our focus, is on the A' couplings. For simplicity and 
feasibility, we keep two non-vanishing A'-couplings. The choice of the mass spectrum and various parameters of the 
model is dictated by the requirement that in the limit of i?-parity conserving SUSY, the branching ratio falls within 
the experimental limits. Then, switching on the non-zero values of the (combination of two) A'-couplings allows 
us to trace their contribution and obtain bounds on the class of RPV parameters from imposing the experimental 
constraints. In this aspect, the strategy is common to many of the earlier studies. We find that RPV contributions 
not only expand the relevant four-quark operator basis of the case of MSSM from 8 to 28, but also introduce new, 
likely to be dominating, contributions in the form of charged-slepton and sneutrino loops to the Wilson coefficients of 
(chromo)magnetic penguins. However, our results clearly show that the new four-quark operators with non- vanishing 
Wilson coefficients at the electroweak/SUSY scale, induce effects through QCD running which typically dominate 
over the direct contributions to the Wilson coefficients of (chromo)magnetic penguins. However, there are exceptions. 
Depending on the particles exchanged, the loop contributions could also dominate in some cases. 

The two-A' type of RPV contributions split into two classes: (A) X'i3jXh2k ^"^^ (-^) ^'ij3^'hk2- Whereas the contri- 
butions in class (A) give rise to hf, —f transition and hence in direct interference with major contributions from SM 
and i?-parity conserving SUSY parts, the contributions in class (B) have not much i?-parity conserving counterparts. 
We obtain numerical bounds, some cases of which are orders of magnitude stronger than available bounds in the 
literature. The interpretation of such bounds has to be taken more carefully. Their values depend strongly on the 
choice of other SUSY parameters. Moreover, it could happen that the i?-parity conserving and RPV contributions 
have some accidental, but strong cancellation. In the latter case, no meaningful bounds on the RPV couplings alone 
can be obtained. After all, the extra contributions in the case of MSSM are expected to have partial cancellations to 
offset the strong charged Higgs contribution. Our numerical results at least indicate clearly the strong implication of 
RPV couplings on 6 ^ s -|- 7 at a level substantially beyond earlier studies are able to illustrate. A most conservative 



24 



interpretation of the bounds obtained would be that they indicate values of the A'-couplings around and above which 
the corresponding RPV contribution alone to & — *■ s + 7 would be of alarming magnitude. We also explain clearly 
why the A"-couplings have no major role to play in & — > s + 7. Another type of significant contributions come from 
combinations of a bilinear and a trilinear RPV couplings. The type of contributions has a quite different structure 
and has not been studied before. They are implicitly included in our comprehensive formulation presented here, while 
the relevant numerical study is to be presented in another publication 24'|. 

All the discussions above, like most of the earlier studies, has not taken into consideration constraints on the RPV 
parameters from the neutrino masses. However, taking the reasonable assumption that all the neutrino masses are in 
the sub-eV range, and that there are no strong cancellations among RPV contributions, very stringent constraints on 
the parameters can be obtained 36. 37,. ,38, ,39). In fact, such constraints are so important that talking about other 
numerical constraints on RPV parameters without reference to the neutrino mass constraints may not be a sensible 
way to take the model seriously. So, we would like to comment on the issue here. Majorana neutrino masses could 
be induced at one loop level by a pair of trilinear RPV couplings A^^^A^^j. Because these violate lepton number 
by two units, they cannot contribute to the decay 6-^5 + 7. However, the case I = m could be used to obtain 
bounds on A^^^. Strong bounds on such individual couplings could significantly strengthen the bounds on some of 
the combinations that contribute to 6 ^ s + 7 . For instance, choosing the spectrum similar to the one for this study, 
one obtains a strong upper bound of jA^ggj — 7.1 x 10~^. The bound on the parameter much less stringent 

(|A^22l = 3-0 ^ 10~^) as it is suppressed by a quadratic factor of {ms/rabf relative to the one-loop contribution from 
the parameter A^gg. This would mean that, whenever the coupling A^„j,„ is a part of the combination contributing to 
& ^ s + 7 , bounds coming from neutrino mass consideration are stronger. Nevertheless, the bounds obtained by us 
on three combinations in case (B), namely |A']^3]^A']^2il7 l'^23i'^22i L I'^331'^32iIj still survive because A^^^ does not figure 
in these combinations. 

note added in proof: Just after this work was finished, an eprint j4Cll | appeared, where the branching fraction for 
B ^ Xs + has been recalculated at the NLL order using the results of soft-collinear factorization for inclusive 
_B-meson decay distributions. It has been pointed out that the singnificant perturbative uncertainties associated 
with the parameter A = mi, — 2Eq have been ignored in the previous works. The new estimate now stands at 
Br[B ^ Xs+-/ {E-y > l.SGeV)]^^ = (3.44 ± 0.53 ± 0.35) x lO"'^, where the first error is the estimate of perturbative 
uncertainties and the second one refiects uncertainties in the input parameters. As a result of this larger theoretical 
uncertainties, the lower bound on the charged-Higgs mass is strongly reduced compared to previous estimates, to 
slightly below 200 GeV at 95 % confidence level. With the larger theoretical uncertainties in the SM prediction, the 
case for new physics certainly becomes stronger. 
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APPENDIX A: THE 7 MATRIX 

In table IIIIHll we provide a translation of our notation of operators to that of refs.poll2l| for easy comparison. We 
define the scheme independent anomalous dimension matrix 7*^^ as: 

7^°^°^ = I " I (Al) 
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where 7°^ and 7^^*^ are given as (for convenience we have given the operators on the top row): 
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TABLE II: Operators that mix in the block 7^ 
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